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While fundamental-mode discrete solitons have been demonstrated with both self-focusing and 
defocusing nonlinearity, high-order-mode localized states in waveguide lattices have been studied 
\^**)" ' thus far only for the self-focusing case. In this paper, the existence and stability regimes of dipole, 

quadrupole and vortex soliton structures in two-dimensional lattices induced with a defocusing 
, nonlinearity are examined by the theoretical and numerical analysis of a generic envelope nonlinear 

■ lattice model. In particular, we find that the stability of such high-order-mode solitons is quite 

different from that with self-focusing nonlinearity. As a simple example, a dipole ("twisted") mode 
i soliton which may be stable in the focusing case becomes unstable in the defocusing regime. Our 

("*\ . results may be relevant to other two-dimensional defocusing periodic nonlinear systems such as 

Bose-Einstein condensates with a positive scattering length trapped in optical lattices. 
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INTRODUCTION 

go : 

Ever since the suggestion of optically induced lattices in photorefractive media such as Strontium Barium Niobate 
(SBN) in 0, and its experimental realization in 0, S , there has been an explosive growth in the area of nonlinear 
waves and solitons in periodic lattices. A stunning array of structures has been predicted and experimentally obtained 
in lattices induced with a self-focusing nonlinearity, including (but not limited to) discrete dipole ||, quadrupole |(|, 
necklace and other multi-pulse patterns (such as e.g. soliton stripes ||), discrete vortices 0, and rotary solitons 
|lfj| . Such structures have a potential to be used as carriers and conduits for data transmission and processing, in the 
context of all-optical schemes. A recent review of this direction can be found in (see also H2)- 
0^ ■ Many of these studies in induced lattices were also triggered by the pioneering work done in fabricated AlGaAs 
waveguide arrays 0|. In the latter setting a multiplicity of phenomena such as discrete diffraction, Peierls barriers, 
diffraction management [lif and gap solitons |l5| among others [lf3 | were experimentally obtained. These phenomena, 
in turn, triggered a tremendous increase also on the theoretical side of the number of studies addressing such effectively 
*sO discrete media; see e.g. nam for a number of relevant reviews. 

Finally, yet another area where such considerations and structures are relevant is that of soft-condensed matter 
physics, where droplets of Bose-Einstein condensates (BECs) may be trapped in an (egg-carton) two-dimensional 
optical lattice potential The latter field has also experienced a huge growth over the past few years, including the 
prediction and manifestation of modulational instabilities [2(| , the observation of gap solitons [2]J and Landau-Zener 
tunneling |22j among many other salient features; reviews of the theoretical and experimental findings in this area 
• »~j ■ have also been recently appeared in [23L |24| . 

In light of all the above activity, it is interesting to note that the only structure that has been experimentally ob- 
served in two-dimensional (2d) lattices in "defocusing" media consists of self-trapped "bright" wave packets (so-called 
"staggered" or gap solitons) excited in the vicinity of the edge of the first Brillouin zone |2J ■ However more complex 
coherent structures have not yet been explored in lattices with defocusing nonlinearity and their stability properties 
have not yet been examined, to the best of our knowledge. It should be mentioned that the defocusing context is 
accessible in the aforementioned settings. E.g., in the photorefractive lattices, this can be done by appropriate reversal 
of the applied voltage to the relevant crystal, while in BECs, the defocusing nonlinearity corresponds to the most 
typical case arising in dilute gases of 87 Rb or 23 Na. 

It is the aim of the present work to examine the non-fundamental soliton structures (e.g., dipoles, multipoles, and 
vortices) in lattices with a defocusing nonlinearity, and to illustrate the similarities and differences in comparison 
to their counterparts in the focusing case. In particular, we study dipole structures (consistingof two peaks) and 
quadrupole structures (featuring four peaks), as well as vortices of topological charge S = 1 (cf. |9j) in a 2D induced 
lattice with a defocusing nonlinearity. These structures will be analyzed in detail for both cases, namely, the "on-site" 
excitation (where the center of the structure is on an empty lattice site between the excited ones) and the "inter-site" 
excitation (where their center is between two lattice sites and no empty lattice site exists between the excited ones). 
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Our study of these structures will be conducted analytically and numerically (in the next two sections) in the 
context of the most prototypical generic envelope lattice model, the so-called discrete nonlinear Schrodinger (DNLS) 
equation with a defocusing nonlinearity which is related to all of the above contexts [l7l l23| . When we find the 
relevant structures to be unstable, we will also briefly address the dynamical evolution of the instability, through 
appropriately crafted numerical experiments. Finally, in the last section, we will summarize our findings and present 
our conclusions, and the interesting experimental manifestations that they suggest. 



MODEL AND THEORETICAL SETUP 



As our generic envelope model encompassing the main features of discrete diffraction and defocusing nonlinearity 
we use the two-dimensional (2D) DNLS equation: 

iii n = -C (A 2 u) n + |w„| 2 u n , (1) 

where w n is a complex amplitude of the electromagnetic wave in nonlinear optics |17| . or the BEC wave function 
at the nodes of a deep 2D optical lattice n is the (two-dimensional in the present study) vector lattice index, 
and A2 the standard discrete Laplacian. Furthermore, C is the constant of the intersite coupling (associated with 
the interwell "tunnelling rate" |23j), and the overdot stands for the derivative with respect to the evolution variable, 
which can be z in optical waveguide arrays, or the time t in the BEC model. We focus on standing-wave solutions of 
the form u n = exp(— iAt)4> n , with </> n satisfying the equation, 

f(cf> n , C) = -H n - CA 2 K + |0 n | 2 n = 0. (2) 

Perturbing around the solutions of Eq. J3J) gives rise to the linearization operator 

W (C) _ f -A + 2|0 n | 2 cbl \ (I 0\ 

with the overbar denoting complex conjugation. Through an appropriate rescaling of the equation, we can fix A = 1. 
Our analysis uses as a starting point the so-called anti-continuum limit, i.e., the case of C = 0, where for the uncoupled 
sites, 

^ = r n e i9 », (4) 

with the amplitude r n being or yA, and the phase 9 a being an arbitrary constant. Continuation of such a solution 
to nontrivial couplings necessitates that a certain, so-called Lyapunov-Schmidt condition be satisfied j2||. The latter 
imposes for the projection of eigenvectors of the kernel of 7ii°' ) onto the system of stationary equations to be vanishing. 
This solvability condition provides a nontrivial constraint at every "excited" (i.e., r n ^ 0) site of the AC limit, namely: 

- 2ig n (6, C) = -Ce- ie "A 2 q> n + Ce w "A 2 ^ n = 0. (5) 

It is interesting (and crucial for stability purposes) to note that this equation has an extra (— ) sign in comparison to its 
focusing counterpart. The derivation of these solvability conditions is especially important because the corresponding 
Jacobian 

Ma = dg./dBj (6) 

has eigenvalues 7 that are directly related to the "regular" eigenvalues of the linearization problem A, through the 
equation 

A = ±v^. (7) 

Hence, the method that we use to derive the eigenvalues A (which fully determine the crucial issue of stability of the 
solution for small C) consists of a perturbative expansion of the solution from the AC limit 

0n = ^ + +..., (8) 

which allows us to derive the principal bifurcation conditions for a specific configuration and therefore infer its linear 
stability properties through the eigenvalues of M. and their connection to the linearization eigenvalues A. Recall that a 
nonzero real part of any eigenvalue is a necessary and sufficient condition for an exponential instability in Hamiltonian 
systems, such as the one considered herein. 
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Type 


On-site 


Inter-site 


Stability 


Instability Outcome 


Stability 


Instability Outcome 


In-phase Dipole 


Unstable 


1-Site Pulse 


C < 0.064 


1-Site Pulse 


Out-of-phase Dipole 


C < 0.092 


Decay 


Unstable 


1-Site Pulse 


In-phase Quadrupole 


Unstable 


Breathing Behavior 


C < 0.047 


1-Site Pulse 


Out-of-phase Quadrupole 


C < 0.08 


1-Site Pulse 


Unstable 


2-Site Mode 


Vortex 


C < 0.095 


1-Site Pulse 


C < 0.095 


1-Site Pulse 



TABLE I: Summary of the stability results for all the configurations presented below. For partially stable (near the anti- 
continuum limit) solutions their interval of stability (for A = 1) is given. In each case, the outcome of the instability evolution 
for the parameters and initial conditions considered below is also mentioned. 

COMPARISON OF ANALYTICAL AND NUMERICAL RESULTS 

General Terminology 

We start with some general terminology that we will use in this section. The designation in-phase (IP) will be used 
for two sites such that their relative phase difference is 0, while out-of-phase (OP) will signify that it is ir. Furthermore, 
on-site (OS) will mean that the center of the configuration is on an empty lattice site (between the excited ones), 
while inter-site (IS) will signify that the center is located between the excited lattice sites (and no empty site exists 
between them). For all modes, in the figures below, we show their power P = Y] \u n \ 2 as a function of the coupling 
strength C, as well as the real and imaginary parts of the key eigenvalues (the ones determining the stability of the 
configuration). We start with the dipole configuration (consisting primarily of two lattice sites; see Figs. 11111) . We 
also examine the more complex quadrupole (see Figs. I5I8[) and vortex (see Figs. 191 10|) configurations. In all the cases, 
we offer typical examples of the mode profiles and stability for select values of C. When the configurations are found 
to be unstable, we also give a typical example of the instability evolution, for a relevant value of the coupling strength. 
Another general feature that applies to all modes is a continuous spectrum band extending for A^ e [A — 8C, A]. This 
latter trait significantly affects the stability intervals of the structures in comparison with their focusing counterparts 
as we will see also below (since configurations may be stable for small C, but not for larger C). 

The presentation of the figures will be uniform throughout the manuscript in that in each pair of figures, we examine 
two types of configurations (one in the left column and one in the right column). The first figure of each pair will 
have five panels showing P as a function of C, the principal real eigenvalues (second panel) and imaginary eigenvalues 
(third panel). In these plots, the numerical results are shown by the solid (blue) line, while the analytical results 
by the dashed (red) line. The fourth and fifth panels show typical examples of the relevant configuration (obtained 
through a fixed point iteration of the Newton type) and its stability eigenvalues (shown through the spectral plane 
(A r , Xt) for the eigenvalues A = A r + «Ai). The accompanying second figure will show the result of a typical evolution 
of an unstable mode, perturbed by a random perturbation of amplitude 10 , in order to accelerate the instability 
evolution. The four contour plot panels (one set on the left and one on the right) will display the solution's squared 
absolute value for four different values of the evolution variable; the bottom panel will show the dynamical evolution of 
the sites chiefly "participating" in the solution. A fourth-order Runge-Kutta scheme has been used for the numerical 
integration results presented herein. 

To facilitate the reader, a summary of the results, encompassing our main findings reported below is offered in 
Table 1. The table summarizes the configurations considered, their linear stability and the outcome of their dynamical 
evolution for appropriate initial conditions in the instability regime. Note that if the solutions are unstable for all C, 
they are denoted as such, while if they are partially stable for a range of coupling strengths, their interval of stability 
is explicitly mentioned. Details of our analytical results and their connection/comparison with the numerical findings 
are offered in the rest of this section. 
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FIG. 1: (Color Online) The first line of panels shows the power P vs. coupling C for the inter-site (IS), in-phase (IP) mode 
(left) and on-site (OS), IP mode (right). The second lines show their maximal real eigenvalues and the third their first few 
imaginary eigenvalues. The solid (blue) lines illustrate the numerical results, while the dashed (red) lines the analytical ones. 
The fourth and fifth panels show the contour plot of the mode profile (fourth panel) and the corresponding spectral plane of 
eigenvalues A = A r + (fifth panel); The left two panels are for the IS-IP mode for C = 0.08 and C = 0.116 respectively. The 
right panel shows the OS-IP mode for C = 0.08. 



Dipole Configurations 

Inter-site, In-Phase Mode 

Figures T1I2I encompass our results for the two types of IP dipole solutions (i.e., initialized at the AC limit with two 
in-phase excited sites) . The IS-IP mode of the left panels is theoretically found to possess 1 imaginary eigenvalue pair 
(and, hence, is stable for small C) 

A « ±2VCi. (9) 

The collision with the continuous spectrum described above causes the mode to become unstable for sufficiently large 
C; the theoretically predicted instability threshold (obtained by equating the eigenvalue of Eq. © with the lower edge 
of the phonon band located at A — 8C) is C — 0.0625, the numerically found one is C « 0.064. Additional instability 
may ensue when the monotonicity of the P vs. C curve changes (we have found this to be a general feature of the 
defocusing branches). The fourth and fifth panels show the mode and its linearization eigenvalues for C = 0.08 and 
C = 0.116. In fact, the dynamical evolution of the mode is demonstrated for the case of C = 0.08, illustrating that 
only one of the two sites eventually persists, after the demonstrably oscillatory instability destroys the configuration 
for t > 100. 
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FIG. 2: (Color Online) The four panels at the top left corner show the contour plots of the dynamical evolution of the unstable 
inter-site, in phase (IS-IP) solution for C = 0.08. The respective times are t = 50 and t = 150 in the top and t — 250 and 
t = 350 in the second row. The panel at the bottom left shows the dynamical evolution of the square modulus of the principal 
two sites participating in the ISTP solution as a function of time. From both of the above, it is clear that the configuration 
relaxes into a single site soliton. The right panels show the same features for the on-site, in-phase (OS-IP) solution, which also 
relaxes (but more slowly) into a single-site configuration. 

On-site, In-Phase Mode 

The OS-IP mode of the right panels of Figs. 1 1121 is always unstable due to a real pair, theoretically found to be 

A rj ±2C, (10) 

for small C . Notice once again the remarkable accuracy of this theoretical prediction, in comparison with the numer- 
ically obtained eigenvalue. The fourth and fifth right panels of Fig. ^show the mode and its stability for C = 0.08. 
Its dynamical evolution in the right column of Fig. shows its slow disintegration into a single-site solitary wave. 

Inter-site, Out-of-phase Mode 

Figures 13141 illustrate the two dipole, out-of-phase modes. The left panels of the figures correspond to the IS-OP 
mode; this one is also immediately unstable (as one departs from the anti-continuum limit), due to a real pair which 
is 

A«2VC, (11) 

for small C. The fourth and fifth panels of Fig. [3] show the relevant mode for C = 0.08 and C = 0.116, showing its 
1 and 2 unstable real eigenvalue pairs respectively. The numerical experiment highlighting the evolution of the mode 
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FIG. 3: (Color Online) Similar to Fig. Q but now for the inter-site, out-of-phase (IS-OP) mode (left panels) and for the on-site, 
out-of-phase mode (OS-OP). The fourth and fifth rows of panels are for C = 0.08 and for C = 0.116 in both cases. 



for the case of C = 0.08 is shown in the left panel of Fig. Clearly, in this case as well, the positive real eigenvalue 
leads to the growth of one of the two sites constituting the dipole, and the eventual formation of a single-site solitary 
pulse. 



On-site, Out-of-phase Mode 

The right panels of Fig. !3l4l show the OS-OP mode. The stability analysis of this waveform shows that it possesses 
an imaginary eigenvalue 

A « 2Ci. (12) 

This leads to an instability upon collision (occurring theoretically for C = 0.1, numerically for C w 0.092) with the 
lower edge (located at A — 8C) of the continuous band of phonon modes. The mode is shown for C = 0.08 and 
C = 0.116 in the right panels of Fig. [3] The direct integration of the unstable solution with C = 0.1 is shown in 
the right panels of Fig. indicating that in this case the mode completely disappears (because of the oscillatory 
instability) into extended wave, small amplitude radiation. 
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FIG. 4: (Color Online) Same as in Fig. H but now for the IS-OP mode (left panels) and the OS-OP mode (right panels). The 
former shown at times t — 25, t = 50 (top row), t — 150 and t = 250 (second row), in the case of C = 0.08 finally results into a 
single-site configuration (as is also illustrated by the bottom panel showing the two principal sites participating in the mode). 
The latter, shown for C = 0.1, at times t — 50, t = 100 (top row), t — 150 and t = 200 (second row) is entirely destroyed by 
the instability resulting into small amplitude excitations. 



Quadrupole Confirgurations 

Inter-site, In-phase Mode 

Figures 15161 show the quadrupolar mode with four in-phase participating sites when centered between lattice sites 
in the left panels of the figures. This mode is theoretically predicted to have two imaginary (for small C) eigenvalue 
pairs with 

A w 1\f<5i (13) 

and one imaginary pair with 

A a V8Ci. (14) 

As a result, this mode (shown in the fourth row panels of Fig. [S]for C — 0.05 and C = 0.1) is unstable due to the 
collision of the above eigenvalues with the continuous spectrum occurring theoretically for C ~ 0.0477, while in the 
numerical computations it happens for C » 0.047. The outcome of the instability shown in Fig. Elfor C — 0.08 is the 
degeneration of the quadrupolar mode into a single-site excitation. 
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FIG. 5: (Color Online) The top three lines of panels show the same features as the corresponding ones of figure Q but now for 
the quadrupole IS-IP mode (left) and the quadrupole OS-IP mode (right). The contour plot of the real part of the modes and 
the spectral plane of their linearization eigenvalues are shown in the fourth and fifth rows for C = 0.05 and C = 0.1 in the case 
of the former mode, while the latter is only shown for C = 0.05. 



On-site, In-phase Mode 



The right panels of Figs. !5l6l show the case of the on-site, in-phase mode. The latter is found to always be unstable 
due to a real eigenvalue pair of 

A « ±4C (15) 

and a double, real eigenvalue pair of 

A ± VUC. (16) 

This can also be clearly observed in the fourth and fifth panels of Fig. showing the mode and its stability for 
C = 0.05. The dynamical evolution of the unstable mode for C = 0.05 is shown in the panels of Fig. Both from the 
contour plots at the different times and from the dynamical evolution of the main sites participating in the structure, 
it can be inferred that the mode embarks in an oscillatory breathing, without being ultimately destroyed in this case. 



Inter-site, Out-of-Phase Mode 

We next consider the case of the IS-OP mode in Figs. 17181 Our analytical results for this mode show that for small 
values of C, we should expect to find it to be immediately unstable due to three real pairs of eigenvalues, namely a 
single one with 

X^±V8C (17) 
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FIG. 6: (Color Online) Same as Fig. but now for the quadrupole IS-IP mode with C = 0.08 (left panels) and the quadrupole 
OS-IP mode with C = 0.05 (right panels). The left panels are for t = 200, t = 300 (top), t = 400 and t = 500 (middle) and 
illustrate alongside the bottom panel (containing the evolution of the principal four sites participating in the structure) how 
the configuration eventually degenerates to a single-site soliton. The right panels are for t — 50, t — 150, (top row) t — 250 and 
t = 350 (second row) and show together with the bottom panel the complex oscillation (breathing) involved in the behavior of 
the quadrupole OS-IP mode for C = 0.05. 

and a double one with 

A ra ±2VC. (18) 

This expectation is once again confirmed by the numerical results of the left panel of Fig. [7| The fourth and fifth 
panels show the mode and the spectral plane of its linearization for the cases of C = 0.08 and C — 0.116. The 
dynamical evolution of this mode also gives an interesting result, in that it produces, upon manifestation of the 
instability, a long-lived, two-site oscillatory mode, as is illustrated in the left panels of Fig. [8] for C = 0.08. 

On-site, Out-of-phase Mode 

Finally, the last one among the quadrupolar modes is the OS-OP mode, examined in the right panels of Figs. 17181 
Our theoretical analysis predicts that this mode should have a double imaginary eigenvalue pair of 

A w ±2Ci (19) 

and a single imaginary pair of 

A re ACi. (20) 

These, in turn, imply that the mode is stable for small C, but becomes destabilized upon collision of the larger one 
among these eigenvalues with the continuous band of phonons. This is numerically found to occur for C ps 0.08, 
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FIG. 7: (Color Online) Similar to Fig. Qbut for the quadrupole IS-OP mode (left panels) and the quadrupole OS-OP mode 
(right panels). The fourth and fifth rows show the modes and their stability for C = 0.08 and C = 0.116 in each case. 



while it is theoretically predicted, based on the above eigenvalue estimates, to take place for C = 0.083. The mode's 
stability analysis is shown in the fourth and fifth panel of Fig. [7| for C = 0.08 and C = 0.116; for C = 0.1, and 
its dynamical evolution is examined in the right panels of Fig. [S] In this case, we do find that the mode essentially 
degenerates to a single-site solitary wave. 



Vortex Configuration 

Inter-site Vortices 

Finally, Figs. 191 101 show similar features, but now for the IS (left panels) and OS (right panels) vortex solutions 
[jfil l2t| . The former has a theoretically predicted double pair of eigenvalues 

A w ±2Ci. (21) 

leading to an instability upon collision with the continuum band for C > 0.095 (C > 0.1 theoretically). In this case, 
there is also an eigenvalue of higher order 

A re ±AC 2 i. (22) 

which obviously depends more weakly on C. The fourth and fifth panels of Fig. [!|]show the real and the imaginary 
parts of the vortex configuration for C = 0.08 and C = 0.116 and the sixth panel the corresponding spectral planes. 
The dynamical evolution of the vortex of topological charge S = 1, for C = 0.1 is shown in the left panels of Fig. ^3 
indicating that the vortex also, upon the occurrence of the oscillatory instability, becomes a single-site solitary wave. 
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FIG. 8: (Color Online) Same as Fig. but for the evolution of the quadrupole IS-OP mode with C = 0.08 (left panels) the 
quadrupole OS-OP mode with C = 0.1 (right panels). The former is shown for t — 25, t = 50 (top row), t = 150 and t = 250 
(second row) substantiating (together with the bottom panel showing the principal four sites of the branch) its resulting into 
a long-lived, breathing two-site waveform. The latter is shown for t = 50, t = 150 (top row), t = 250 and t = 350 (second row) 
indicating its degeneration into a single-site configuration. 



On-site Vortices 

On the other hand, the OS vortices are shown in the right panels of Figs. 191 101 In this case, we theoretically find 
that the vortex, for small C, should have a double pair of eigenvalues 

A w 2Ci (23) 

and a single, higher order pair of eigenvalues 

A re ±V32C 3 i. (24) 

The former eigenvalue pairs, upon collision with the continuous spectrum, lead to an instability, theoretically predicted 
to occur at C — 0.1 and numerically found to happen for C pa 0.095. The on-site mode (and its stability) is shown 
in the fourth-sixth right panels of Fig. [5] for C = 0.08 and C — 0.116. Its evolution (for C — 0.1) is shown in the 
right panel of Fig. ^3 where it is again seen that the mode degenerates from an S = 1 to an S = structure, i.e., a 
single-site solitary wave with no topological charge. 



General Principles Derived From Stability Considerations 



It is interesting to note as an over-arching conclusion that the stability intervals of the defocusing structures are 
different from those of their focusing counterparts (especially when they are stable close to the AC limit) because 
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FIG. 9: (Color Online) The same features as in Figure 1 are shown here for the IS vortex of topological charge S = 1 (left) 
and the OS vortex of S = 1 (right). In this case both the real (fourth line) and imaginary (fifth line) parts of the solution are 
shown (and their stability in the sixth line) for C = 0.08 and C = 0.116. 



of the collisions with the continuous spectrum band edge; the latter is at A = A in the focusing case, while it is 
at A = A — 8C in the defocusing setting. Another similarly general note is an immediate inference on whether the 
structures are stable or not; this can be made based on the knowledge of whether their focusing counterparts are stable 
or not and the transformation from the former to the latter through the staggering transform: u n , m = (— l) n+m v n , m - 
For instance, IP two-site configurations (both OS and IS) are known to be gencrically unstable in the focusing regime 
[2(|; through the staggering transformation, OS-IP of the focusing case remains OS-IP in the defocusing, while IS- 
IP of the focusing becomes IS-OP in the defocusing. Hence, these two should be expected to be always unstable, 
while the remaining two (OS-OP in both focusing and defocusing and IS-OP of the focusing, which becomes IS-IP 
in the defocusing) should similarly be expected to be linearly stable close to the AC-limit, as is indeed observed. 
Notice that, interestingly enough, for the vortex states the staggering transformation indicates that the stability is 
not modified between the focusing and defocusing cases. This is because for an IS vortex, it transforms an S = 1 
state into an S = — 1 state (which is equivalent to the former, in terms of stability properties), while the OS vortex 
remains unchanged by the transformation. However, as mentioned above, these considerations are not sufficient to 
compute the instability thresholds for initially stable modes, among other things. They do, nonetheless, provide a 
guiding principle for inferring the near- AC limit stability of the defocusing staggered states, based on their focusing 
counterparts. 



CONCLUSIONS AND FUTURE CHALLENGES 



In this paper, we have studied in detail some of the principal multi-site solitary wave structures that emerge in 
the context of defocusing nonlinearities, examining, in particular, dipole, quadrupole and vortex configurations. We 
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FIG. 



t t 

10: (Color Online) Same as in Fig. |21 but for the dynamical evolution of the IS (left panels) and the OS (right panels) 



vortex of topological charge S = 1. Both cases are shown for C = 0.1. Both the left and right panels of evolution are for t = 50, 
t — 150 (top row), t = 250 and t = 350 (second row) and alongside the bottom panels showing the principal four sites of the 
vortex they show how the vortices transform themselves into fundamental solitary waves centered on a single site (in the latter 
case, interestingly, the site of the largest excitation of the ensuing wave is not one among the four principal excitation sites of 
the original OS vortex). 



have illustrated which ones among these states can potentially be stable (e.g. IS-IP and OS-OP for both dipoles 
and quadrupoles, as well as the vortices) and those that will always be unstable (e.g. IS-OP and OS-IP modes 
for both dipoles and quadrupoles). We have also provided detailed analytical estimates of the stability eigenvalues 
associated with these modes, in very good agreement with the observed numerical results. The analytical calculations 
also empower us to identify, even for the stable (close to the AC-limit) modes, the relevant intervals of stability 
of those waveforms. We have corroborated our analytical calculations with detailed computations that identify the 
corresponding modes and numerically analyze their linear stability. In addition, for each of the modes, we have shown 
some typical examples of their dynamical evolution, when they become unstable (either directly, or subsequently due 
to eigenvalue collisions). 

These results offer immediate suggestions for experiments in arrays of optical waveguides, Bose-Einstein condensates 
(e.g. of 87 Rb or 23 Na, which feature repulsive interactions amounting, at the mean-field level, to a defocusing 
nonlinearity) mounted on a deep optical lattice. In the latter case, the nodes of the lattice considered herein would 
correspond to BEC droplets in the respective wells of the optical potential. Finally, they are also suggestive of 
similar experiments in the recently and rapidly growing theme of photorefractive crystal lattices (where, however, the 
nonlinearity is slightly different, featuring a saturable form). 

We close by suggesting that these results also indicate that higher charge configurations pjl . \F\ may similarly be 
possible and could potentially also be stable in a defocusing setting, similarly to the 5 = 1 states discussed above. It 
would certrainly be of interest to examine such states in the near future, as well as to study the effect of additional 
components |28| (i.e., multi-component states, relevant to the above optical settings when multiple polarizations are 
present, or to BECs when multiple hyperfine states are studied), or that of higher-dimensional structures [29j- 
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